Non-analytic spin susceptibility of a nested Fermi liquid: the case of Fe-based 

pnictides 
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We propose an explanation of the peculiar linear temperature dependence of the uniform spin 
susceptibility x(T) in ferropnictides. We argue that the linear in T term appears due to non-analytic 
temperature dependence of x{T) m a two-dimensional Fermi liquid. We show that the prefactor of 
the T term is expressed via the square of the spin-density-wave (SDW) amplitude connecting nested 
hole and electron pockets. Due to an incipient SDW instability, this amplitude is large, which, along 
with a small value of the Fermi energy, makes the T dependence of x{T) strong. We demonstrate 
that this mechanism is in quantitative agreement with the experiment. 

PACS numbers: 71.10.Ay, 75.30.Cr, 74.25.Ha, 74.25.Jb 



Introduction. The key hypothesis of the Fermi liq- 
uid (FL) theory is that a system of strongly interact- 
ing fermions can be considered effectively as a gas of 
weakly interacting quasiparticlcs. In the absence of resid- 
ual interaction between quasiparticles, the static uniform 
spin susceptibility x{T), and the specific heat coefficient, 
7(T), are finite at T = and obey quadratic dependen- 
cies on T at low temperatures. The effect of residual 
interactions on x(T) and j(T) has been studied inten- 
sively in recent years [l[ , with the key result that in two 
dimensions (2D) both j(T) and x{T) are linear rather 
than quadratic in T @, y, 

Theory predicts that the behavior of 7(T) is univer- 
sal in a sense that the (negative) slope is given by the 
square of the backscattering amplitude [Bj|. This linear 
decrease has been observed in monolayers of 3 He @. On 
the contrary, a linear in T term in the spin susceptibility 
is not universal and can be of either sign 0, 0] , causing 
uncertainty in the interpretation of the experiments on 
semiconductor heterostructures [1]. 

Recently, a pronounced linear temperature dependence 
of the uniform susceptibility has been observed in high- 
T c su per conductors with iron-based layered structure 
Hell 111 EH • It extends from temperatures above either 
SDW or superconducting transitions up to 500-700K with 
almost doping-independent slope. The T dependence is 
quite strong: x(T) increases roughly by a factor of 2 be- 
tween 200K and 700K. An explanation of this behavior 
based on the Ji — Ji model of localized spins has been 



The itinerant character of Fe-pnictides is suggested by 
the agreement between the band structure obtained in 
ab initio calculations [lit ano - observed in de Haas-van 
Alphen and angle-resolved photoemission spectroscopy 
(ARPES) experiments 0, G3, GjJ H3 . It is firmly 
established by now that the Fermi surface (FS) of Fe- 
pnictides consists of two small hole pockets near (0, 0) 
and two electron pockets near (tt, ir) points of the folded 
Brillouin zone. In such a system, an obvious origin of 
the SDW order is a logarithmic divergence of a particle- 
hole vertex involving states on nested parts of the FS 
HmSHIll Si- 



proposed in Ref. [13|; however, given that the linear T de- 
pendence persists up to large dopings, where local probes, 
such as nuclear magnetic resonance (NMR) and /iSR, do 
not see localized moments 14, [Hj], this explanation is 
questionable. 



In this Letter, we propose an explanation of the ex- 
perimental T-dependence of spin susceptibility based on 
the itinerant picture. We argue that the origin of the 
linear increase of x{T) with temperature in ferropnic- 
tides is the same as in a 2D FL. Furthermore, we show 
that this behavior is universal for FLs with strong (n, ir) 
SDW fluctuations, namely, the slope of the linear in In- 
dependence is determined by the square of the SDW am- 
plitude with nesting momentum Q = (7r,7r). This am- 
plitude is large which, along with a small value of the 
Fermi energy ep, amplifies the T dependence of xC^)- 
Choosing the SDW coupling to reproduce the observed 
SDW transition temperature Tjv at zero doping, wc find a 
good agreement between calculated and measured slopes 
of x(T). The doping dependence of the slope also agrees 
with the data. We view this agreement as a strong in- 
dication that the linear temperature dependence of x(^) 
in pnictides @, 03, 11, T^] is in fact the first unambigu- 
ous observation of a non-analytic behavior of the 2D spin 
susceptibility. Besides being fundamentally important on 
its own right, this observation also strengthens the case 
for the itinerant scenario for Fe-pnictides. 
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FIG. 1: Second-order diagrams for the thermodynamic poten- 
tial. Solid and dashed lines correspond to /-fermions (elec- 
trons) and to c- fermions (holes), respectively. 



Theory. We consider a two-band model of interacting 
fermions occupying the electron and hole FSs: 
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Here, Ck<j (/ker) is the annihilation operator for a hole 
(electron) with momentum k and spin a (for an electron, 
k is measured from the (ir, 7r)-point), e k and e£ = — £ k + 
represent single-particle dispersions, and measures a 
deviation from perfect nesting. Models of this type were 
considered jn the past in the context of an "excitonic 
insulator" 



27] 



We assume that each of the electron and hole FSs is 
doubly degenerate. The terms with U4 and 1*5 are intra- 
band interactions, the terms with u\ and 112 are inter- 
band interactions with momentum transfer and Q, re- 
spectively, and the term with U3 is the inter-band pair 
hopping. All couplings flow from their initial values at 
energies of order of the bandwidth to rcnormalizcd values 

" 5F 



at £f 
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We assume that this rcnormalization is 



already included into Eq. (fTJ) and analyze the behavior of 
the system at energies below ep, where the actual values 
of momenta become relevant. 

We first obtain the linear-in-T contribution to the spin 
susceptibility, &x{T), to second order in the interaction 
and then show that the prefactor of the T term in Sx(T) 
is expressed via the SDW vertex to all orders in the in- 
teraction. 



Fig. [T] depicts all topologically incquivalcnt second- 
order diagrams for the thermodynamic potential, 
$(T, H), describing both intra- and inter-band processes. 
Each of the diagrams contains the object 
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where G£ u is a Matsubara Green's function of either c or 
/ fcrmion with momentum k and spin-cr in a magnetic 
field H . As in an ordinary FL, the non-analytic H 2 T 
term in <p(T,H) comes from a dynamic Kohn anomaly, 
i.e. from diagrams with momentum 2kp carried by the 
interaction lines 0, 13] ■ It is, however, more conve- 
nient to re-express the result via the polarization bub- 
bles with small rather than 2kp momenta. The Green 
functions in Eq. ^ can be combined in two different 
ways: either as n^IT^', or as n^n^, where 11^ = 

T E^I (0 G k,c G £+q^+o is a polarization bubble of 
the cc, //, or cf type with small momenta q <C hp- A 
bubble depends on the magnetic field if the Zeeman en- 
ergies of two fermions add up. One can readily show that 
the cc and // bubbles depend on the field via n , while 

the field enters the cf bubble through n^ n and 
terms. Evaluating individual diagrams, we find that each 
of them can be expressed as a product of two dynamic 
spin up/down bubbles: 



Hq,n = 



2^ ^{n-2ifi B H -2i A i*) 2 + {v F q) 2 ' 



(3) 



where /1* = for intra-band scattering and /1* = \i 
for inter-band scattering. The rest of the calculations 
proceeds in the same way as for an ordinary FL 0, |3] • 
In short, one integrates II q n over q first, replaces the 
Matsubara sum by a contour integral, differentiates the 
result twice with respect to H, sets H = 0, and ob- 
tains the 0(T) terms in %. The first two diagrams in 
Fig. [1] describe intra-band processes and give the same 
results for X (T) as in Refs. d fl: b~XiA T ) = Xu l.5 T i 
where A = 4xo-Nf,/(2e F ), e F = vpkp/2, \o = 2fi%N F 
is the Pauli susceptibility per one sheet of the FS, and 
Np = m/2n is the density of states. A factor of 4 in 
A results from the double degeneracy of electron and 
hole bands. The third through fifth diagrams involve 
inter-band scattering and yield Sx3{T) = 2\u 2 Tr)(n/T), 
5 Xi (T) = X^Tnin/T), and 6 X5 (T) = 2Xu 2 T V (p/T), re- 
spectively. Here, tj(x) = 2xcothx — x 2 /sinh x — 2x. The 
third and forth diagrams give finite contributions if a = 
a', while the fifth diagram contributes if a = —a' . The 
sixth and seventh diagrams do not contribute to S X (T), 
and the remaining two give 8xs{T) = —2XuiU2Tri(^/T) 
and 8X9(7) = -Xu 2 T V (v/T). 



Combining all diagrams, we obtain 

S X (T) = XT [u\ + u\ + 2 (u\ + u\ - uiu 2 ) vi^/T)] . 

(4) 

The first two terms are the contributions from intra- 
band processes - they are the same as in an ordinary 
FL. The rest of the terms correspond to inter-band pro- 
cesses, specific to the electronic structure of nested fer- 
ropnictidcs. For T 3> U (perfect nesting), r](fj,/T) 1 
and both intra- and inter-band processes contribute to 
the linear term in the spin susceptibility. In the opposite 
limit of T < fx (poor nesting), r?(/i/T) ~ exp(-2u/T) 
and the inter-band contribution to x(T) is suppressed 
exponentially. In the intermediate regime T ~ fi, rj(fi/T) 
is non-monotonic. Note that the pair-hopping term M3, 
which gives rise to an attraction in the extended s-wave 
(s±) pairing channel, does not contribute to the T term 
in X (T). 

The coupling constants in Eq. ((4| are interactions at 
the scale of the Fermi energy. These couplings are al- 
ready renormalizcd from their bare values at the scale 
of the bandwidth [24| by parquet renormalization group 
(RG), in such a way that U4, M5, and U\ flow to the same 
value, while u 2 flows to a smaller value. The renormal- 
ized coupling at the scale of ep depend only weakly on 
the incoming and transferred momenta. Further renor- 
malization below ep differentiate between Ui with differ- 
ent transformed and incoming momentum. Such renor- 
malization is particularly relevant for our system as the 
coupling ui, which corresponds to the scattering process, 
Mi (k, p + q; p, k + q) 4/^ +q / p+q Ck, diverges at the on- 
set of the SDW instability for q = (we recall that mo- 
menta for /-fermions are measured from Q, so that the 
momentum transfer between c and / fermions is actu- 
ally q = Q). The singular vertex is iti(k, p; p, k), which 
means that electrons and holes swap their respective mo- 
menta. Note that the divergence occurs for any angle 
between k and p. At weak coupling, the enhancement of 
Mi is confined to q <C kp, while for a generic q ~ fcp the 
coupling Mi retains its bare value. 

We now show that the u\ term in x(T) oc T is the 
same coupling that diverges at the SDW instability. To 
see this, we first note that the fermionic momenta in 
diagrams for &(T,H) are constrained by two require- 
ments: i) the momentum transfers are near 2kp and ii) 
all four momenta are near the FS. For the third diagram 
in Fig. [TJ this implies that p ss — k, |k| sa kp, while q is 
small (~ T/vp). Therefore, the vertex in this diagram is 
Mi(k, — k; — k, k). This is an analog of the backscattering 
amplitude in a 2D FL. The vertex Mi(k, — k; — k, k) is a 
special case of the SDW vertex Mi(k, p; p, k) for p = k. 

Next, we consider higher-order diagrams. They can 
be separated into two classes. In diagrams of the first 
class, one obtains the non-analyticity by keeping only 
two dynamic bubbles n q o and lumping the rest of the 
diagram into renormalization of the static scattering am- 
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FIG. 2: (a) T dependence of the spin susceptibility x(T), as 
given by Eq. (|6]), for perfect nesting (n — 0) and for a range 
of couplings uNf, as indicated in the plot. In a calculation, 
max(T, /i) in Eq. (|6]) is approximated by y^T 2 + ^t 2 . Inset: 
Calculated slope dx(T)/dT as a function of uN F at T = 300K 
(chosen to match Ref.Qj]). (b) Calculated x(T) for uN F = 0.5 
for a range of /1, as shown in the plot. Inset: Same curves as 
in (a) normalized by xifyextr, obtained by extrapolating y(T) 
down to T = 0. (c) x{T) in BaFe 2 - a; Co ;c As 2 from Ref. [ij. 
Inset: same data as in the main panel normalized by x(0)- (d) 
Calculated slope dx(T — 300K) jdT as function ofu. Inset: 
measured slope as a function of doping (from Ref. \lw . 



plitudc. In particular, the backscattering amplitude 
Mi(k, — k; — k, k) is renormalized into an effective cou- 
pling uf 1 , which diverges at the SDW instability. Sin- 
gular renormalizations of Mi form ladder series which is 
summed into 



.off 



Ml 



1 - uiNpln 



(5) 



max(T,/i) 



The second class of higher-order diagrams contain three 
and more dynamic bubbles. Terms of order of m^ etc. are 
not expressed in terms of backscattering amplitude but 
rather contain Mi with typical q of order kp. In this range 
of momenta, Mi is not enhanced by SDW fluctuations and 
remains small at weak coupling. 

Neglecting U2 and setting ui = M4 = M5 = u > 0, we 
finally obtain 



Sx(T) 



8 {uNpf xo 
vpkp 



T 



V(J*/T) 



1 - uN F In 



max (T, p;) 



(6) 

We see that the full result for 6 X (T) is obtained from the 
second-order expression by replacing mi = u by the exact 
SDW amplitude Mf ff , given by Eq. . This is similar to 
the result for the specific heat [5j, but differs from the 
result for 6x(T) in an ordinary 2D FL, where the full 
Sx(T) is not expressed via the backscattering amplitude 
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|3|, [Z| . This difference can be traced down to the symme- 
try between the particle-particle (Cooper) channel in an 
ordinary FL and the particle-hole channel in a nested FL. 
Indeed, the backscattering contribution to both j(T) and 
Sx(T) for the ordinary case undergoes logarithmic renor- 
malization in the Cooper channel. However, the Cooper 
ladder for the ordinary case is identical to the particle- 
hole ladder for the nested case, except for that the sign 
of the interaction is reversed, i.e. the SDW instability 
for u > for the nested case is related to the Cooper 
instability for u < for the ordinary case. In both cases, 
there are also non-backscattering contributions 5\{T). 
For the ordinary case, the backscattering term is reduced 
by Cooper renormalization, and non-backscattering con- 
tributions play the dominant role. For the nested case, 
SDW renormalization enhances the backscattering term, 
and other contributions can be ignored. On the other 
hand, the Cooper channel composed of electrons and 
holes for the nested case is equivalent to the particle- 
hole channel for an ordinary case and, therefore, is not 
logarithmically divergent. 

Comparison with experiments. We now apply Eq. ([6]) 
to ferropnictides. The experimental results for x(^) m 
BaFe2_ 2: Co a ;As2 [12] are shown in panel (c) of Fig. (2] 
From the data, we estimate the slope of the T dependence 
as (x(700K)/x(0)) cxp ~ 2, where x(0) is obtained by ex- 
trapolating x{T) to T = (theoretically, x(0) ~ 4xo)- 
Taking v F = 0.45eV-l and k F sa 0.16A" 1 from the 
ARPES data [IH, we obtain e F ~ 0.04 eV 0. The 
only unknown parameter of the theory - the dimcn- 
sionless coupling constant uNp - is fixed by requiring 
that the SDW vertex uf ff increases upon approaching 
T N {T N = 140K at zero doping). As Eq. (JSJ is an ap- 
proximate one- loop formula, we set a criterium that u e f f 
increases by a factor of 2 at Tjy. This yields uNp w 0.5. 
We then find ( X (T = 700K)/x(0)) thcOT « 1.7, which is 
quite close to the experimental (x(700K)/x(0)) exp ~ 2. 
A more detailed comparison between the experiments 
and Eq. (|6|) is presented in Fig. [2l where we also show 
the dependencies of the slope on uNp and fi. We find 
quite a good agreement with the experimental data. 

Conclusion. To summarize, we analyzed a non- 
analytic, linear in T term in the spin susceptibility of 
a 2D Fermi liquid with nested electron and hole pockets 
of the Fermi surface. We found that the prefactor of the 
T-term contains the same inter-band coupling , which 
is enhanced by SDW fluctuations. These results describe 
quantitatively a strong temperature dependence of the 
spin susceptibility in ferropnictides, observed in a num- 
ber of recent experiments. An immediate consequence of 
the proposed mechanism is that x should exhibit equally 
strong linear dependencies on the magnetic field and on 
the wave number Q . We suggest to perform these mea- 
surements as a crucial test for the origin of the observed 
effect. 
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